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Abstract Nomenclature 


The Rotor Systems Research Aircraft 
uses load cells to isolate the 
rotor/transmission system from the fuse- 
lage. An analytical model of the rela- 
tionship between applied rotor loads and 
the resulting load cell measurements is 
derived by applying a f orce-and-moment 
balance to the isolated rotor /transmi ssion 
system. The model is then used to esti- 
mate the applied loads from measured load 
cell data, as obtained from a ground-based 
shake test. Using nominal design values 
for the parameters, the estimation errors, 
for the case of lateral forcing, were 
shown to be on the order of the sensor 
measurement noise in all but the roll 
axis. An unmodeled external load appears 
to be the source of the error in this 
axis. 


Introduction 


The Rotor Systems Research Aircraft 
(RSRA) has a set of seven load cells con- 
necting the main rotor transmission to tha 
fuselage. Their purpose is to make high- 
accuracy measurements of the net rotor 
loads, as resolved at the rotor hub, from 
flight data (Ref. 1). The use of these 
load cells to estimate applied rotor 
forces and moments at the hub requires an 
accurate mathematical expression relating 
rotor loads, inertial loads, and load cell 
readings. Both the structure and param- 
eters of this model must be specified. 
Previous approaches to processing ground- 
test data have not yielded an acceptably 
accurate relationship. This is particu- 
larly true for the case of applied high 
frequency dynamic loads. This paper 
describes a new approach to obtaining the 
relationship, and presents the results of 
a preliminary evaluation of the resulting 
model from experimental data. 


A, B,C,D, E, F,G - individual load c'll 
outputs ( Fig . A1 ) 

a * linear acceleration at 

C.G. 


b - position of C.G. from 

shaft attach point 

d - longitudinal distance 

between vertical load cell 
attach points 

E{ } - expected value operator 


e - lateral offset of for 

ward lateral load cell 
from centerline 


e H 
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H 

ft 

h 

I 


error in estimate o “ 
applied hub loads 

lateral offset of aft 
lateral load cell from 
centerline 

vector of six applied 
hub loads 

estimates of applied hub 
loads 

position of hub from C.G. 

moment of inertia about 
C.G. 


J 


vector of six inertial 
loads 


1 - position of C.G. from 

aft load cell attach 
points 

M -massirat.rix 


m - effective mass of 

rotor/cransmission/engine 

ays ten 

N 9 - number of samples 
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Subscripts 

A 


= index of measurement 
samples 

« rotational rates about 
x , y , z axes 

■ covariance of measure- 
ment noise 

- total applied engine and 
tail rotor drive torques 

- transformation matrix 
for hub loads 

» transformation matrix 
for load cells 

- vector of seven load 
cell measurements 

- lateral distance between 
vertical load cell attach 
points 

* force and moment compo- 
nents (Fig. A2 ) 

* gyroscopic coupling 
coef f i ci en t s 

- set, of unknown parameters 

* a priori estimate of 
parameters 

- post-calibration esti- 
mate of parameters 

- measurement noise 

« load cell deformation 
angles about x f y, and 
z axes, respectively 

- frequency 


« accelerometer 


a ,b , c , d , e , f ,g * attach point of each 
load cell 


Previous Methods 

The initial attempt to determine the 
load cell response to the applied rotor 
load3 ‘nvolved applying static loads at 
the hub and measuring the resulting load 
cell response. A least-squar os regression 
approach was used to identify a coeffi- 
cient matrix relating the seven load cells 
to the six hub loads. A six element bias 
vector was also estimated. It was found 
(Ref. 2) that the coefficient matrix 
varied as a function of the applied mul- 
tiple axis load. This indicates a non- 
linear dependence of the load cell 
response to the applied rotor loads, and 
would require a polynomial expansion of 
the multi-input multi-output relationship 
to characterize it in terms of constant 
paramet ers . 

The next calibration was a ground- 
based shake test in which a pair of iner- 
tial shakers were mounted on the rotor hub 
to apply dynamic loads at specified ampli- 
tudes and frequencies. During this exper- 
iment, the RSRA was suspended from the hub 
so the static loading was the same for all 
tests. As a result, the nonlinear varia- 
tion of the relationship with static loads 
observed in the static test should not be 
present in the shake test. The Force 
Determination Method (Ref. 3) was utilized 
to estimate applied rotor loads from a 
variety of sensors around the aircraft. 
This method first identifies the transfer 
functions between the applied rotor loads 
and sensors at various points on the air- 
craft, then identifies tne applied rotor 
load from a least . quares fit to the 
transfer functions and measured sensor 
responses. The results were unacceptable 
because the identified transfer functions 
varied with the magnitude of the applied 
load; hence, they could not be used to 
estimate the applied load witnout a more 
extensive calibration procedure. Since 
there was only one static load condition, 
it appears that this nonlinearity is due 
to a different mechanism than the nonline- 
arity observed in the static tests. 


C - total load 

H - hub load 

I -inertial load 

m • measured data 

T - load cell 

x,y,z -component for x,y,z axis 


Pro p osed Metho 

Since the sensors utilized in the 
Torce Determination Method (FDM) included 
numerous accelerometers and strain gauges 
mounted on the fuselage, transfer func- 
tions of these sensors will be affected by 
any nonlinear dynamic behavior in the 
fuselage. This effect complicates the use 
of fuselage sensors to determine applied 
rotor loads. The RSRA was designed to use 
load cells to isolate applied loads from 
diffrerent sources, such as the main 
rotor, tall rotor, engine, and wings. The 
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proposed approach takes advantage of this 
concept by treating the rotor transmission 
as an isolated system with externally 
applied loads from the load cells and the 
rotor (Fig. 1). The applied rotor loads 
are then measured from a f or ce-and-moment 
balance using measured load cell loads (T) 
and inertial loads (J) as derived from 
transmission acceleration measurements. 

In order to utilize this approach, a model 
is required chat relates applied external 
loads on the rotor /transmi ssion system to 
the resulting forces and moments at the 
center of mass of the system. This model 
is derived analytically from jhysical 
principles, using the known eometry of 
the rotor /transmissi on system. Parameters 
with potentially uncertain values in this 
model are explicitly represented to pro- 
vide the capability to calibrate the 
modei . 

The advantages of this approach pre- 
dominantly arise from the physical insight 
obtained in using an analytically derived 
model. With such a model, sensitivity 
analysis and physical judgment can be used 
to select the m o y t appropriate set of 
available parameters for calibration. 
Consequently, fewer parameters need be 
calibrated than when no physical insight 
is used. In addition, the parameters to 
be calibrated now have a physical inter- 
pretation so that the validity of the 
calibration results may be assessed. The 
model should initially be derived to be as 
simple as possible. If it cannot ade- 
quately explain the observed experimental 
behavior with a physically reasonable set 
of parameter values, it may be expanded to 
include additional effects, as required. 

It is important that all major effects be 
identified and Incorporated into the model 
before calibration of the parameters is 
attempted, or the parameter values wi'll 
compensate for the unmodeled effects as 
best they can and achieve physically unre- 
alistic values in the process. 


Objective and Approach 

The objective of this study is to 
derive a simple dynamic model of the iso- 
lated rotor transmission system and test 
its accuracy with experimental data. 

The approach is first to derive a 
simple model of the rotor transmission 
system, trsating it as a linear, rigid, 
isolated body. Known or assumed values 
are used for a^l parameters of the derived 
model. The model is then applied to test 
data to determine its accuraoy. It it 
appears that calibration can further 
improve the accuracy, the appropriate 
^a'ameter set will be selected and cali- 
brated. If the accuracy appears limited 


by an unmodeled effect, the model will be 
expanded to evaluate potential sources of 
the unmodeled effect. 


Rotor /Transmission Model 

The arrangement of the load cells 
below the rotor/transmission system is 
shown in Fig. 2. A detailed description 
of this system, including all inertia 
contributions, is given in Ref. 2. There 
arc seven load cells; four are mounted 
vertically at the corners of the trans- 
mission mounting plate, two are mounted 
laterally at the fore and aft edges of the 
mounting plate, and one is mounted longi- 
tudinally at the forward edge of the 
mounting plate. An inertial ' oad vector, 
J, is located at the center * cavity of 
the rotor/transmission system, the 

applied rotor load vector (H) is located 
at the rotor hub (Fig. 3). The rotor hub 
is located at the end of the rotor shaft 
which is tilted forward at an angle of 
2°. The load cells are connected to the 
transmission and the fuselage by spherical 
bearings . 

fhe proposed approach is to estimate 
the applied hub loads from a force-and- 
moment balance of the external loads and 
the inertial loads. In order to accom- 
plish this, all externally applied loads 
must be transformed to the center of mass 
where the inertial loads act. By treating 
the rotor /transmi ssion as an isolated 
system, the load cell forces are consid- 
ered a measured, externally applied load 
on the system. A 6 x 7 matrix (S) is 
derived that transforms the seven load 
cell loads (T) at their attach points to a 
set of six load components at the center 
of mass. A 6 x 6 matrix (R) is derived 
that transforms the applied rotor load (H) 
at the hub to the center of mass. An 
inertial load (J) at the center of mass is 
derived from measured accelerations and 
assumed inertial parameters. Using the 
derived matrices, a f or ce-and-moment 
balance at the center of mass results in a 
set of Si v simultaneous equations which 
may be wri -.ten in matrix form as: 

J+S*T'R*H»0 (1) 

A detailed description of these vectors 
and matrices are given in Appendix A. 

The assumptions used in deriving the 
matrices and Eq. (1) are that the 
rotor /transmi ssion s>stem is a rigid body 
and that there is no friction in the load 
oell bearings. These assumptions we*e 
made to simplify the Initial approach. 

Both nonrigid body effects and friction in 
the load cell bearings could be added to 
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the model in order to match the experi- 
mental data if that appears to be 
required . 

The transf or mat 1 on matrices S and 
R and the inerti il load vector J were 
all derived to explicitly contain all 
potentially uncertain parameters of the 
system, so that any subset of parameters 
can be selected for calibration. The 
parameters defined in the model are: 

a) All distances and angles required 
to define the resultant moment arm from 
the center of mass to load application 
points . 

b) All angles and magnitudes 
required to define load cell load compo- 
nents acting on the rotor /transmission 
system . 

c) All mass properties required to 
determine inertial loads from measured 
rates and accelerations. 


Test Cond i t i ons 

Having derived a model, the next step 
is to evaluate it with experimental 
data. Both static and dynamic ground test 
data are available. The dynamic data 
generated by the shake test were chosen 
since they would provide a more rigorous 
test of the model structure than would the 
static test data. The dynamic data are 
not, however, suitable for testing the 
calibration procedure. This is because 
the same static load condition exists for 
all dynamic tests, and parameter varia- 
tions are mostly dependent on variations 
in the static loading. Once the model 
structure has been validated with the 
dynamic data, the model can be applied to 
the static test data to evaluate the cali- 
bration procedure. 

The test datum selected was a fre- 
quency sweep from 15 to 18.5 Hz in the 
y-axis applied rotor load (lateral 
force). This frequency range was chosen 
because it contains the N/rev frequency, 
and identification of applied rotor loads 
at this frequency is of special 
interest. The lateral forcing was chosen 
because previous tests have shown the 
poorest results with y-axis forcing, so 
it would provide the most rigorous test. 
Transfer function data were generated from 
the raw test data by a harmonic analyzer 
for four levels of applied load. Equa- 
tion (1) was then used to generate the 
applied load estimate from the transfer 
function data. Equation (1) was processed 
with all of the model parameters set to 
assumed nominal values. 


Error Analysis 
T otal Estimation Error 

In a controlled ground test environ- 
ment, the actual /alues of the applied 
rotor loads are known so the total error 
in the estimate can be readily obtained. 
Given the measured load cell readings, Tm , 
and inertial loads derived from accel- 
erometer measurements, Jm, Eq . (1) may be 
used to estimate the applied rotor loads 
as : 

H - -R 1 * [Jm + S • Tm] (2) 

The total error in the estimate is 
obtained by subtracting the known values 
of applied rotor load from the estimate of 
Eq . ( 2 ) to get : 

e^ e H - H * -R • [Jm + S * Tm] - H 

(3) 

e H - -R' 1 * [Jm + S • Tm ♦ R • H] 

The available shake test data were in the 
form of transfer functions that had been 
generated from the raw data by a harmonic 
analyzer. In order to utilize these data 
with the proposed model, it was necessary 
to transform the model to the frequency 
domain and write it in terms of the trans- 
fer functions. Transforming Eq , (3) to 
the frequency domain gives: 

e H U) 

- -R ^ • [Jm[u>) + S • Tm(w) + R * H(w)] 

m 

If only a y-axis rotor load, Hy(w), is 
applied, Eq. (4) may be written in terms 
of the transfer functions as: 

- -R 1 • [ ( Jm( w) /Hy ( u>) ) 

+ P * (TmUj/HyU)) + R] • Hy(w) 

(5) 

where Jm(to)/Hy(oJ and Tm ( w) /Hy ( <u) are 
vectors of transfer functions of the iner- 
tia) and load cell loads with respect tc 
the y-axis rotor load. Equation (5) is 
us-id to evaluate the total error In the 
estimate using the available transfer 
function data. 
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Numerous potential sources of error 
are present in this system. The evalua- 
tion procedure is to examine the total 
error and attempt to categorize it into 
the potential sources. Once identified, 
the sources would be modeled and included 
in the system. Most of the error will 
probably be attributable to one of five 
sources : 


A) Systematic errors in data collec- 
tion . 

B) Unmodeled static and dynamic 
effects . 

C) Unmodeled external loads. 

D) Random errors in sensors. 

E) Incorrect parameter values. 
Systematic Errors in Data Collection 

One source of error is the use of 
transfer function data. Since this is 
treated as raw data in this study, any 
errors in the identification of the trans- 
fer functions would propagate through the 
proposed approach. The recorded time 
domain data should be reprocessed by the 
harmonic analyzer to provide only Fourier 
transformed data, not transfer functions. 

Angular acccl erometer s and rate gyros 
were not available on the rotor transmis- 
sion system for the shake test. Since 
there is no way to obtain all such data, 
the approach taken is to assume it is 
negligible and see how the estimates com- 
pare with this assumption. Some justifi- 
cation for this assumption comes from 
comparing the response of the two linear 
accelerometers mounted on the transmission 
with a 2-ft vertical displacement between 
them. The difference in the y-axis com- 
ponents divided by the vertical displace- 
ment should give the roll axis angular 
acceleration. The average value of this 
dei j.ved roll acceleration over the fre- 
quency range for an applied load of 800 lb 
was found to be 0.09 rad/sec/sec, support- 
ing the low angular acceleration assump- 
tion. The derived angular acceleration 
data was not used with the model because 
the errors in the linear accelerometers 
are such that the accuracy of the derived 
angular acceleration is 0.6 rad/sec/sec. 
The derived values are therefore in the 
noise level. The angular acceleration 
affects the translational equations since 
the linear accelerometer is not mounted at 
the center of mass and will, therefore, be 
affected by angular accelerations. The 
moment equations are affected since the 
angular inertial loads are dependent on 
the angular accelerations. The errors 
associated with neglecting a 


0.09 rad/sec/sec roll angular acceleration 
are 12 lb in the y-force and 27 ft-lb in 
the rolling moment so the assumption 
appears justified. 

Unmodeled Static and .Dynamic Effects 

The model derived for this study was 
kept deliberately sample to facilitate the 
analysis. It can be expanded if necessary 
to account for the observed error. Unmod- 
eled effects with potentially significant 
Impact on the model include friction in 
the load cell bearings, flexibility in the 
rotor/transmi ss ion system and nonlinearity 
in the dynamic response to applied 
loads. The nature of the error signal 
should suggest which c: these effects are 
>resent. Friction and deadbands will be 
characterized by hysteresis in the 
response. This effect is more readily 
observed in static data than in dynamic 
data. Flexibility will show up as a reso- 
nance at some frequency and will result in 
a phase and amplitude shift between the 
input and output signals. Nonlinear 
dynamics will be readily detected by a 
frequency shift between the input and 
output data. Static nonlinearities result 
in parameter variations and are corrected 
by calibration rather than by expanding 
the model. 

Random Errors In Sensors 


Both load cells and accelerometers 
have measurement noise that produces a 
lower bound on the accuracy of the applied 
load estimates. It is possible to obtain 
accuracies below this limit, but this 
requires the use of statistical processing 
techniques such as Kalman Filtering and 
Smoothing. This effect can certainly not 
be reduced by any modifications to the 
model. The effect of accelerometer and 
load cell noise on the applied rotor load 
estimate is derived in Appendix B and used 
to generate the errors given in Table I. 
These numbers were based on the assumption 
of independent random errors for each 
sensor with accuracies of 1 % of full scale 
for the accelerometers and 0.1J of full 
scale for the load cells. 

Unmoidled External Loads 


The derived model will be in error if 
all externally applied loads are not 
included. If the levels are low, L-hen 
this can be the most difficult source of 
error to identify. This is because it can 
take on virtually any ohar ac ter i s t i c* 1 and 
will blend in with other rrror sources. 

The only possibility for detecting this 
type of error is if it is sufficiently 
large that it cannot be logically 
explained by any of the other error 
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sources. Once the presence of an unmod- 
eled mad is suspected, the error data may 
assist in isolating its source, but a 
thorough examination of the test condi- 
tions is usually required to resolve this 
effect . 

Incorrect Parameter Values 

once the model structure has been 
validated to the fullest extent possible, 
the remaining errors should be due only to 
incorrect values for the parameters. At 
this point, calibration may be applied to 
reduce this error source. If calibration 
is attempted before the model structure is 
adequately determined, the parameters will 
take on whatever values are required to 
compensate for the model structure 
errors. This will result in physically 
unrealistic values for the parameters and 
could, in fact, be a test for whether the 
model structure is accurate. 

The conventional approach to calibra- 
tion has been to apply least-squares mini- 
mization of the error with respect to the 
parai.eters to be calibrated. Since the 
parameters are now imbedded in a model, 
the least-squares minimization of the 
error must be done subject to the con- 
straint that the model equations are 
satisfied. This is refered to as a con- 
strained least-squares approach (Ref. *1) 
and the algorithm is derived in Appendix C 
for the constraint of Eq . (1). 

Since the available data are in the 
frequency domain, the calibration must be 
performed in the frequency domain. This 
is actually an advantage since the fre- 
quency domain transformation has concen- 
trated information for the required fre- 
quency range into fewer data points than 
required for a time domain representa- 
tion. Calibration may then be performed 
witn fewer data : n the frequency domain. 
The leas t -squares minimization nay be 
applied to frequency domain data in the 
same way as it Is applied to time domain 
data (Ref. 5). The only modification is 
that the data are organized with the real 
and imaginary parts stacked end to end 
rather than using the data in complex 
form. 'I' his insures that the identified 
parameter values will not be complex. 


Results 

Equation (3) was applied to transfer 
function data for the load cells and the 
main rotor gear cox accelerometers to 
generate the estimation error of the 
derived model in the frequency domain. 
Assumed, nominal \alues, based on aircraft 
design specifications, were used for all 
parameters in the model. Figures H and 5 


show the magnitude of the transfer func- 
tion data for the seven load cells and two 
accelerometers for an 800 lb y-axis exci- 
tation. The location of the load cells 
may be seen by finding the correspondingly 
labeled load cell in Fig. 2. The acceler- 
ometer transfer functions have been multi- 
plied by the s item mass co produce iner- 
tial loads. The strong coupling in the 
system is apparent from the large vertical 
load cell values for a y-axis excita- 
tion. The strong correlation in the 
z-axis is particularly apparent from the 
sudden drop in z-axis acceleration at the 
same frequency (16,6 Hz) where two verti- 
cal load cells (A and B) suddenly assume 
equal and opposite values. 

Figures 6 through 11 show the error 
i.i the model, as defined by Eq. (3), using 
nominal parameter values. Table I shows 
the average error over the frequency range 
as compared to the accuracy limit set by 
the instrumentation noise. With the 
exception of applied rolling moment, the 
average error shown in Table I and the 
frequency plots shown in Figs. 6-11 demon- 
strate that the rigid body model of the 
isolated rotor transmission system pro- 
duces applied load estimates with error 
levels comparable to the instrumentation 
noise level. 

The pronounced roll moment error 
(Fig. 9) is too great to be explained by 
parameter errors, sensor errors, or angu- 
lar accelerations. The demonstrated lack 
of significant angular accelerations about 
the roll axis Indicates that it is not due 
to a nonrigid body effect. It can also be 
seen that there is no frequency shift (the 
hump at 16.6 Hz in the error signal 
matches the hump in the load cell data at 
that frequency) so the error <s not a 
nonlinear function of the modeled vari- 
ables. The only remaining explanation is 
that the error is due to an unmodeled 
extern: 1 load with predominant effect in 
the roll axis. A potential candidate for 
the source of this load is the drive train 
sinoe it would affect only the roll 
moment. Measurements of the shaft torques 
from ti.e engine and the tail rotor are 
required to verify this and could be used 
to ctmpens&te for this effect. 

Nominal parameter values produce 
estimation errors on the order of the 
instrumentation noise level in all but the 
roll axis, and the roll axis error is too 
great to be explained by parameter uncer- 
tainties! therefore, calibration is not 
needed to improve the accuracy in five 
axes i d would not nelp in the sixth. Th 1 
sou"c* of unknown roll axis loading mus* 
be d*:c*rmined and modeled before the r ; l 

>e used for roll moment estimation. 


38U 





Conclusions 


Appendix A: Load Cell Model Derivation 


An analytically derived model with 
nominal parameter values has been used to 
estimate applied rotor loads from measured 
load cell and accelerometer data. This 
approach has also provided a check on the 
consistency of measured input/output 
data. The presence of an unmodeled exter- 
nal load was detected in the roll axis and 
its source is being Investigated. 

The following conclusions are drawn 
from this work: 

For the y-axis case, an analytically 
derived linear, rigid body model of the 
isolated rotor/transmissi on system with 
nominal parameter values performs well in 
all but the rol 1 axis. 

Nonrigid body effects or nonlinear 
behavior canno f explain the roll moment 
error. The error must be due to an exter- 
nally applied, unmodeled load. 

Calibration is unnecessary f r this 
case since improved parameter estimates 
will not improve the estimation error 
further. The estimation error in all but 
the roll axis is already in the noise 
level and the error in the roll axis is 
too great to be logically explained by 
different parameter values. 


A diagram of the undeformed l<~ad cell 
geometry is shown in Fig. A1 . The load 
cells are mounted to the fuselage and „o 
the transmission base by spherical bear- 
ings. It is assumed that the transmission 
base dees not warp, so all changes in the 
load cell geometry are due to deformations 
in the load cells themselves or in the 
fuselage mounting points. With this 
assumption it is possible to completely 
model the load cell response using a gen- 
eral three-component representation of the 
reaction force at each attach point on the 
transmission base; no knowledge of the 
fuselage deformation is required. 


Three plane views of the load cell 
geometry with the t hr ee -component reaction 
force representation are shown in Fig. A2 . 
Inertial loads, including gravity and the 
drive shaft torque, are assumed to be 
concentrated at the transmission center Oi 
gravity. The rotor loads are concentrated 
at the hub. Taking the sum of the forces 
and moments about the center of gravity 
gi ves : 


X T + X + G + E + F + A 
I H x x x x 

+ B + C + D - 0 
xxx 


Recommendations for Further Research 
The following recommendations are 

made : 


Y + Y 
I H 


G + ? + F * A 

y y y y 


♦ B + C ♦ D - 0 

y y y 


Apply the model to data generate by 2 
excitation in other axes to validate It 
under a broader range of conditions and to 
further isolate the source of the unmod- 
eled load. This analysis should be backed 
up by a follow-on physical error source 
analysis. 

L C 

Apply the model to static test data 
to evaluat ■ the calibration technique. 

Nominal parameter values have done sur- 
prisingly well in allowing accurate hub 
load estimation for the specific static 
load condition present in the dynamic 
data. 5i tatic test experience suggests 
that this will not be true under a 1*1 
static load conditions. Once the model 
structure is fully validat'd it shouLd be M C 
applied to static teat data to determine 
which parameters to calibrate and to 
determine the range of variation of the 
parameters through calibration. 


Z+Z+G + E + F + A 
I H z z z z 

♦ B + C + D - 0 
z z z 


L i + l h * y h 


h ♦ ( A + D ) 
z z z 


w/2 


E z • e ' {B Z + C z ) 
-(B + C + F + G 

y y y y 


w/2 + F 


E + A ♦ D ) 

y y y 


m x + m h 


h - Z„ • h 

z H x 


- (B Z + A z + h * V •<«-■»> 


+ (C x * D x + F x + B x + \. 
. b t ♦ (C z * D z ♦ F z ) • 1 


E + G ) 
x x 
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°y * V 
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1 ) - (C + F + D ) 

y y y 


(A x * V 


+ V 


w/2 - E 
w/2 - F 


x 

f + Y 


e + (C 


Using the transformations: 


- A cos (♦ a )sln( 0 a ) 
A y = -A si n ( ) 

A^ - A cos ( <> a ) cos ( 0 a ) 
= B cos f <j> b ) si n ( 0 b ) 


8^ = -B si.i(4» b ) 

B - B cos ( ) cos ( 0. ) 

z Y b b 

C * C cos C 4> ) s in ( C t ) 

X CO 

C « - C s I n ( ) 

y c 


C - C cos ( )cos ( 0 ) 
z c c 

D * D cos ( <J> . ) 3 i n ( 0 . ) 
x d d 

Dy - -D sin(<^ d ) 

D « D cos ( <J> . ) cos ( 0 . ) 
z d d 


E x « -E cos ( ) s i n ( 4> e ) 

Ey = E cos ( $ e ) cos ( i|> e ) 

E » E si n ( ) 

Z 6 

- F cos ( ) s i n ( ) 


Fy * -F cos ( )cos ( 4> f ) 
F z * -F ain($ f ) 

G x - G cos ( 0^ ) cos ( i|/ g ) 

G - G ccs { 0 ) s i n ( * ) 

0 - -G a i n ( 8 ) “ 

z g 


the equations become 

J + R*H>w » T * 0 


0 

( A 1 ) 


CA2) 


wher e : 


and 


[X I’ Y l’ Z I’ L r M I* N I ] 


( AM) 


(A3) 
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X T 


-m • a 

I 


X X 

y t 


-m • a + r • q 

I 


y y y 

Z T 


-m • a 

I 


z z 

4 


-Ix • p + Qt + r L • q 

4 


-Jy • q ♦ r M • p 

4 


- Iz * r 



- 


(A5) 


In Eq. (A5), a x , a , and a z arc linear 
accelerations in each direction measured 
at the ro t or / transmi ss 1 on system center of 
gravity, and p, q, and r are rotational 
accelerations. Because of the nonrigid 
engine mountings, the engine contributions 
to inertial forces are not equal in all 
directions when measured at the system 
center of gravity. This effect can be 
adequately modeled by assigning different 
values to the total effective mass of the 
combined rotor/transmission/engine system 


in each direction: 


m x , m y , and m x . There 
are also a few minor error terms riot given 
here that are discussed fullv in Ref. 2. 
Ix, Iy, and Iz are moments of inertia; 
related terms in the cross-products of 
rotations* . ates are negligible. Qt is 
total applied engine and tail rotor shaft 
torque. Gyroscopic coupling forces due to 
engine and transmission lOtational moments 
of momentum are represented by the coeffi- 
cients r, with subscripts fcr the appro- 
priate axis. The hub forces (H) and load 
cell readings (T) in Eq . (A3) are given 
by : 

H* ■ [Xjj , Y^, Mp| , N^] (A6) 

T ' - [A, B, C, D, E, F , G] (A7) 


The geometric transf ormationa 
applied loads to the center of 
S ) are given hy ; 


1 

0 

0 

0 

-h. 


1 0 
0 1 

4 0 

0 -h v 


* om 1 1. e 
i.ass (R and 


( A8) 


V-i 

* 


to 


386 



ORIGINAL PAGE IS 
OF POOR QUALfTY 


and 


C 0 aSea 
- S 0 a 


C 0 bS 0 b 

C0CS0C 

C0dS9d 

-C 0 eS 0 e 

“C 0 f S 0 f 

c*gc*g 

-S 0 b 

-S0C 

-S 0 d 

C0eC'>3 

C 0 fC 0 f 

C>gs*g 


S 


C 0 aC 0 a 


C 0 bC 0 b 


C 0 cC 0 c C 0 dC 0 d 


S 0 e 


30 f "30g 


w/2C0aC0a 

+bzS 0 a 


bzS 0 b 

-w/ 2 C 0 bc 0 b 


bzS 0 c 

- w/ 2C 0 cC 0c 


b z S 0 d 

+w / 2C0dC0d 


-eS 0 e 

-bzC 0 eC 0 e 


-fS 0 f 
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+bzC 0 f S0f 
+ 1S0f 
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+ ( d- 1 ) S 0 g 

(d-1 )C 0 gS 0 g 


(A9) 


where C and S denote cosine and sine 
functions, respectively 


Appendix B: The Effect of Sens o r Noise 

on Load Estimation 


V H - -R _1 • lM 


v A + S • v T ] l B7 ) 


The analytical model of the 
rot^r/transmission system derived in the 
text has the form: 

J+S * T + R • H - 0 (Bl) 


Substituting Eq. (PI) into Lq . (B 6 ) then 
gives: 

H - H + v H , v> H - M(0,Q h ) (B 8 ) 


Given the measured load cell -ead- 
ings, Tm, and inertial loads derived from 
accelerometer measurements, Jm, the 
applied rotor loads, H, are estimated from 
Eq. (Bl) as: 

H - -R" 1 • [Jm + S • Tin] (82) 

The measured load oell and aocelerometer 
values may be written in terms of their 
true values and a random measurement noise 
component as follows. 


am - a + , 

v A - N( 0 ,Q a ) 

(B3) 

Tm - T ♦ v T , 

v T - N(0 ,Ci t ) 

( 3" ) 


Writing the derived inertial load vector, 
Jm, as the product of an inertia matrix, 

M, and the acceleromet er measurements, am, 
gi ves 

Jm - M * am * M • a ♦ M • v, - J + M 

A A 

(B5) 

where J is the actual inertial load 
vec;o . 

Substituting Fqs. ( (B5) into 
Eq. ( B2 ) gives: 


H • -R -1 • [J ♦ S • T] ♦ v H <B6) 


where 


From Eq. (B 8 ) it is seen that v„ , ac 
given by Eq. (B7). is the combined effect 
of the instrumentation errors on the 
applied rotor load estimate, H. This 
error reprise. ts a lower bouna on the 
accuracy of the estimation that is attain- 
able without applying statistical process- 
ing, suoh as Kalman filtering or smooth- 
ing. The covariance of this error, Q H , 
may be oomputed from Eq. (B7) using the 
known covariance of the Instrumentation 
errors, Q A and Q T , as follows. 


Q 

E{ v 

H^H> 




■ 

' E{[ 

•r" 1 

• (M • v A 

♦ r 

v T )] 


• c 


• (M • v A 

♦ s 

• v T ) 3 1 ) 






(B9) 
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i VT” 


Vl 


j'-i 

Vi 




- 1 


The sensor noise components ar^ assumed 
independent so the term iS 

zero. Substituting the known covariance 
matrices of the sensors for the other 
expected value terms gives: 


Q» - R 


-1 


Cm q a . m 

A 


♦ S 


s’ ] • (R' 1 ) 


( B 1 1 ) 


Appendix C: Calibration Algorithm 


The applied rotor load estimation 
error is gWen by Eq . (3) in the text as: 


-H 


• [ Jra + S 


Tm 


H] (Cl) 


The coefficient matrices, R and S 
contain geometric parameters of the model 
and the derived inertial loads vector, Jm, 
include inertial parameters of the 
model. A subset, © , of the parameters is 
selected for calibration and the error is 
treated as a function of those param- 
eters. A cost function is written in the 
f orm : 


N 


vc 0 ) 


1 / 2 


l e ( 0 » n ) . W 

n-1 H 


-1 


■e H ( e » n ) 


( C2 ) 


where W is a weighting matrix given by: 


W - 1 /N • E e„( fl,n) -e u ( 0,n) (C3) 

9 n-1 H H 


Expanding the first partial of V with 
respect to the parameter set 0 about 
their ncninal values of *0 gives: 


3 V( 6) / 3 0 . 3 v( 0) / 3 0 + (O - 0) • 3 2 V(0)/30 2 


(Ci() 


Setting the desired value at 0 to zero, 
corresponding to an extremum, gives; 


6.0- ( 3 2 V( 0) / 3 Q 2 ) “ 1 • (3V(0)/30) 


( C5 ) 


Taking the first partial of V with 
respect to 0 gives; 


av( 6)/30 - l e (6,n) 
n-1 


W • 3e ( ©, n ) / 3© 
H 


( C6 ) 


The second partial, to first order, is 
then : 


N 


3 2 V(0)/30 2 - l ( 3e ( 0,n ) / 30) 1 • W 
n-1 H 


-1 


•(3e H (0,n)/3e) 


(C7) 


Once the parameter set has been chose* » 
the partial derivative of the errors, e^, 
with respect to the parameters must be 
generated analytically from Eq . (Cl). the 
parameter estimates are then obtained from 
Eqs . (C5) to ( C 7 ) using the measurement 

error sequence, e H , generated by Eq . (Cl). 
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Fig. 5 Magnitude of inertial load 

transfer functions for 800 lb. 
applied lateral force. 
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Fig. 7 Magnitude of y-axis applied 
force estimation error. 
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Fig. 10 Magnitude of y-axis applied 
moment estimation error. 
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DISCUSSION 
Paper No. 24 


EVALUATION OF A LOAD CELL MODEL FOR DYNAMIC CALIBRATION OF THE ROTOR SYSTEM RESEARCH AIRCRAFT 

R. W. Du Val 
M. Bah rani 
and 

B. Wellman 


Charlie Fredrickson, Sikorsky Aircraft : I’d like to thank you for your paper and for mentioning 

the drive train and its possible impact on rotor loading, transmission loading, etc. One effect 
that you may not have modeled seems to be kind of glaring and looking at the analytical model in 

your paper was the engine. I see the engine kind of cantilevered off the transmission [or] the 

drive shaft and I am sure that Is not the way it is supported. If it does have a longitudinal 
restraint to the airframe that may account for part of the modeling error that you may not have 
accounted for. 

Du Val : Yes, in fact the way the engine is mounted is that there are swivel mounts in a univer- 
sal Joint that should prevent any longitudinal restraint on the engine. Now Just how well those 

are working is another question. 

Dev Banerlee, Hughes Helicopters : Ron, a very logical and systematic modeling approach to 

identifying the model of the transmission mount. When you get to correlating the test results 
with your analytical model, I presume you would include the noise characteristics into your 
error function. How do you Intend putting in the covariance of the noise of the load cell? 

[How do you] intend to find it? 

Du Val : That basically was done analytically. Since we have a covariance for the load cells 

and for the accelerometers we can use our analytical model to determine what the resulting 
covariance would be in the estimation of the hub loads, using those assumed covariances for the 
sensor data, and in fact that is the way we arrived at the noise level lines that we put on the 
charts . 



